Modelling the post-Newtonian test-mass gravitational wave flux function for compact 

binary systems using Chebyshev polynomials. 
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We introduce a new method for modelling the gravitational wave flux function of a test-mass par- 
ticle inspiralling into an intermediate mass Schwarzschild black hole which is based on Chebyshev 
polynomials of the first kind. It is believed that these Intermediate Mass Ratio Inspiral events 
(IMRI) are expected to be seen in both the ground and space based detectors. Starting with the 
post-Newtonian expansion from Black Hole Perturbation Theory, we introduce a new Chebyshev 
approximation to the flux function, which due to a process called Chebyshev economization gives 
a model with faster convergence than either post-Newtonian or Pade based methods. As well as 
having excellent convergence properties, these polynomials are also very closely related to the elu- 
sive minimax polynomial. We flnd that at the last stable orbit, the error between the Chebyshev 
approximation and a numerically calculated flux is reduced, < 1.8%, at all orders of approxima- 
tion. We also find that the templates constructed using the Chebyshev approximation give better 
fitting factors, in general > 0.99, and smaller errors, < 1/10%, in the estimation of the Chirp mass 
U ' when compared to a fiducial exact waveform, constructed using the numerical flux and the exact 

^ : expression for the orbital energy function, again at all orders of approximation. We also show that 

• in the intermediate test-mass case, the new Chebyshev template is superior to both PN and Pade 
approximant templates, especially at lower orders of approximation. 

> ■ 

VTi ' I. INTRODUCTION 

^ ! 

, The inspiral of compact objects such as neutron stars and black holes are expected to be a major source of 
(3 ■ gravitational waves (GW) for the ground-based detectors LIGO, VIRGO, GEO600 and TAMA [lllllli, as well as 
, for the future planned space-based detector LISA Q . Due to the effect of radiation reaction, the orbit of the binary 
■ system slowly decays over time. As this happens the amplitude and frequency of the waveform increases emitting a 
' 'chirp' waveform. 

There have been many efforts to create templates which will approximate a possible signal to high accuracy. 
^ On one hand we have the post-Newtonian (PN) expansion of Einstein's equations to treat the dynamics of the 
(^JT) system This works well in the adiabatic or slow- motion approximation for all mass ranges. On 

L * the other hand we have black hole perturbation theory El El IH El IH IH which works for any velocity, but only 
• ^ [ in situations where the mass of one body is much greater than the other. While templates have been generated to 5.5 

• PN order for a test-mass orbiting a Schwarzschild black hole E9i ^-nd to 3.5 PN order for non-spinning binaries of 
H I comparable mass ^2^3; ^ number of difficulties still need to be tackled. The main problem is that both templates 

are a function of the orbital energy and GW flux functions. In the test-mass case, an exact expression in known for 
the orbital energy, but we have a PN expansion for the flux function. In the comparable-mass case, a PN expansion 
is known for both functions. It has been shown that the convergence of both methods is too slow to be useful in 
creating templates that can be confidently used in a GW search 0, 0, |2l l2ll . We also know that the PN 
approximation begins to break down when the orbital separation of the two bodies is r < lOM |22j|. This means that 
as we approach the Last Stable Orbit (LSO) the templates begin to go out of phase with a possible signal due to 
the increase of relativistic effects. As most search methods are based on matched filtering, any mismatch in phase 
between our templates and a signal will result in a loss of recovered signal-to-noise ratio (SNR) and an increase in 
the error in the estimation of parameters. 

It was shown |2l |3 |2l| that templates based on resummation methods such as Pade approximation have a faster 
convergence in modelling the gravitational waveform. The Pade based templates were then used to partially construct 
Effective One Body templates |2l which went beyond the adiabatic approximation and modelled the waveform 
into the merger phase. Other more phenomenological templates such as the BCV (2l |2l [s^ templates seem to be 
excellent at detecting GW, but are not necessarily the best template to use in the extraction of parameters. 

In this paper we focus on the inspiral of IMRI sources. These sources encompass the inspiral of a Neutron star (NS) 
into a black hole with masses of 10s to 100s of solar masses (which should be observable in the ground based detectors), 
to the inspiral of a low mass Supermassive black hole (10^ — 10** M0) into a more massive black hole (10^ — IO^Mq) as 
should be observable with future space based detectors. We don't believe that the PN approximation used here will 
be sufficiently accurate to model EMRI sources, and expect that other methods such as analytic and numerical kludge 
waveforms |^^,|3l0| will be used to properly model EMRI waveforms. On the other hand, we fully believe that 
this method of resumming the PN series using Chebyshev polynomials will be also applicable to the comparable mass 



2 



case [35|- 

Throughout the paper we use the units c = G = 1 . 

A. Improving Template Construction Using the Minimax and Chebyshev Polynomials. 

The problem with expansions Uke a Taylor series is that they are based on Weierstrass's theorem, which assumes 
that there are enough terms in the expansion to sufficiently model the function we are approximating. We know from 
previous studies that the 11 term expansion for the flux function for test-mass systems may not be sufficient. A more 
promising possibility is based on getting close to the minimax polynomial by using the family of Ultraspherical (or 
Gegenbauer) polynomials which are defined by 

= C„ (1 - ^ (1 - (-1 < a < ex.) , (1) 

where C„ is a constant. These polynomials are orthogonal over x S [—1,1] with respect to the weight function 

(l — x^)" . A feature of the polynomials P^\x) is that they have n distinct and real zeros and exhibit an oscillatory 
behaviour in the interval [—1,1]. For a = —1/2 the amplitude of the oscillations remain constant throughout the 
interval and is conducive to trying to find an "equal-ripple" error curve, which is integral to the minimax polynomial. 

This value of a corresponds to the Chebyshev polynomials of the first kind, Tn{x), (hereafter Chebyshev polyno- 
mials). These polynomials are closely related to the minimax polynomial due to the fact that there are n + 1 points 
in [-1,1] where Tn{x) attains a maximum absolute value with alternating signs, i.e. |T'„(a;)| = ±1 [s^. It can be 
shown jSTj l that the Chebyshev polynomials exhibit the fastest convergence properties of all of the Ultraspherical 
polynomials. 

For our purposes, we need to approximate polynomials which are a function of the dimensionless velocity v in the 
domain v e [Q.viso], where viso — l/-\/6 is the velocity at the LSO for a test-particle orbiting a Schwarzschild black 
hole. In this case we use the Shifted Chebyshev polynomials, designated T*{v). We can also transform from the 
interval [1,-1] to an arbitrary interval [a, 6] using 



V a; e [a,6],s G [-1,1]. (2) 



2x- {a + b) 

S = 7 

— a 

In this case we have 

s = — -l = V24f-l \/ V e[0,viso]. (3) 

Vlso 

We can now write the shifted Chebyshev polynomials in the form 

t:{v) = r„(s) = r„(V24« - 1), (4) 

and the recurrence relation as 

T:{v) = 2{V2Av - 1)T:_,{v) - T:_,{v), (5) 
such that the shifted polynomials have the initial conditions 

T*iv) = l, T*{v) = V24v-l. (6) 

II. THE TEST-MASS GRAVITATIONAL WAVEFORM. 

In the stationary phase approximation the Fourier transform for positive frequencies is given by [ssl Is^ l40l l4ll] 

/oo 
h{t) exp(27^^/^) dt = ^/-V6g4^(/)-f ] ^ (7) 
-oo 

where A is a, normalization constant. The phase of the Fourier transform in the stationary phase approximation, 
ip{f), is found by solving the set of 1** order ODEs given by 

dt TTTO^ E'if) 

df df 3vj F{f) 
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where m is the total mass, Vf = (ttto/)^/'^ is the instantaneous velocity, E'{v) — dE/dv is the derivative of the orbital 
energy with respect to the velocity v — (mfi)^/"^ — x^^^, where is the angular velocity as observed at infinity and x 
is an invariant velocity parameter observed at infinity. Finally, Fiv) is the gravitational wave flux function. 

For a test-mass particle in circular equatorial orbit about a Schwarzschild black hole, an exact expression for the 
orbital energy exists |42|. Its derivative with respect to the velocity is given by 



E'{v) — —rjv 



1 - 



(l_3„2)3/2' 



(9) 



where we have introduced a finite-mass dependence through the reduced mass ratio, rj = mim2/m?. From this 
equation the LSO is found by demanding E'(v) = 0, giving viso — For the flux function we only have a FN 

expansion of the form 
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where Fiq {x) is the dominant Newtonian flux function given by 

F^{x) = f ^^1°, 
5 



(10) 



(11) 



and the coefRcients in the expansion of the flux function are given by 0, 0, ITHl IT^ UtI IT^ . We begin to encounter 
logarithmic terms at fc = 6 and above. It is well know that terms such as these can destroy the convergence of a 
power series expansion. 

The Fade approximant to the flux is defined as 



FpM = 1- 



'^pole 



F, 



N 



1 + ln 



V 
Vlso 
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fc=6 
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E h-' 



.k=0 



(12) 



where P^j is the Fade operator, Vpoie = l/"\/3 is the velocity at the photon ring and the coefficients Ik and fk are 
related to the coefficients in the original FN expansion. 



III. THE INTERMEDIATE TEST-MASS CHEBYSHEV APPROXIMATION FLUX FUNCTION TO 

5.5-PN ORDER. 



After factoring out the logarithmic term and introduce a linear term into the non-logarithmic series, the next step 
is to expand both power series in the above equation in terms of the shifted Chebyshev polynomials. This is done 
by writing each monomial in the power series in terms of the shifted Chebyshev polynomials and substituting back 
into the series above. So, for example, starting with Equations ((SJ and 0, we can invert each expression for the 
monomials in w, i.e. 



l = T*{v), v^{2A)-^/'[T;{v)+TI{v)] 



(13) 



and so on. 

Froceeding like this for all monomials, it then allows us to write the power series in the FN expansion solely in 
terms of the shifted Chebyshev polynomials. The first advantage the Chebyshev approximation has over both the 
FN and Fade approximations is that we can also expand the power series appearing in the logarithmic terms as a 
Chebyshev series. Substituting for the monomials, we can write 
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E = E ^^T^niv), fy = kt:{v). (m) 

fc=6 fe=0 

As the expression for both series in terms of the coefficients , A^. are long so we will omit them here. We should 
emphasise the fact here that although the values of the coefficients Ik are zero up to fc = 6, the Chebyshev expansion 
includes terms from k = 0. This allows us to define the Chebyshev approximation to the gravitational wave flux 
function as 



FcM = 1 



'^pole 



Fn{v) 



Vis 
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1 + ln — >.4T:W 



k=0 
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E 

J Lk=0 



(15) 



where we re-introduce the pole at the photon ring. 
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n 


5 


6 


7 


8 


9 


10 


11 


^C„ 


4.3 X 10"^ 
9.9 X 10"* 


4.7 X 10"^ 

1.8 X 10"'' 


7.8 X 10"^ 
7.6 X 10"^ 


5.4 X 10"^ 
2.6 X 10"^ 


1.7 X 10"^ 
5 X lO"'' 


1.4 X 10"^ 

8.5 X 10"^ 


2.8 X 10"^ 
7.1 X 10"* 


tot 
Tn 
tot 

c„ 


1.37 

1.3 X 10"^ 


0.94 
2.9 X 10"" 


0.47 
1.1 X 10"'' 


0.39 
3.2 X 10"^ 


0.34 
5.9 X lO"** 


0.17 

9.2 X 10"^ 


2.8 X 10-^ 
7.1 X 10"* 



TABLE L The top two lines give the truncation errors associated with each order of approximation for both the PN and 
Chebyshev flux functions at the LSO, where the error in the approximation is greatest. The bottom two lines give the total 
truncation error incurred as we reduce the order of approximation from 5.5 to 2 PN. 



A. Modelling the test-mass flux function. 

A major reason for using Chebyshev polynomials is a process called Chebyshev economization )43l] . This works as 
follows : we can always expand a power series Pn{x) of some function f{x) in terms of a Chebyshev series qn{x) such 
that 

n n 

Pn(x) = ^ flfcx'' = ^ hkTk{x) = qn{x). (16) 

A:=0 fc=0 

In general the series Pn{x) diverges as we move away from the point of expansion. Also, we usually find that 
\oin\ > |art-i|. Therefore if we wish to truncate the polynomial pn{x) to give another polynomial rm{x), such that 
m < n, we introduce an additional truncation error of the order |a„a;"| in rm{x). On the other hand, the coefficients 
hk in a Chebyshev expansion always decrease as we increase the order of approximation such that |6ri| < |&n-i|- The 
same truncation with the Chebyshev series qn(x) to give another Chebyshev series Sm(x), gives an error which is 
guaranteed to be smaller than an equivalent truncation in the Taylor series. This is due to the fact that as |T„(a;)| < 1 
in the interval, the truncation error is of the order of the truncated coefficient |6„| which is in general a small number. 
The other main advantage of using Chebyshev polynomials is the fact that because they use information from the 
entire interval, their convergence has absolutely nothing to do with the convergence of the initial Taylor series. 

It has been shown in the past that for a test-mass in a circular equatorial orbit around a Schwarzschild black hole, 
that resumming the PN flux approximation using Fade approximation produces a better flt to the numerical flux. 
However, this method does suffer from certain aspects such as poles and singularities, and the fact that it may not 
perform as well as the PN flux at low orders of approximation. In this section we investigate the modelling of the 
numerical flux |4^ . calculated using the Teukolsky formulism from black hole perturbation |45l |. using Chebyshev 
economization. It has also been well known that the PN approximation has its largest error at the LSO. We would 
expect as v takes on a maximum value at the LSO, that the maximum truncation error should be well approximated 
by the first neglected term, i.e. for Ft„_i{x) the error should be smaller than 

eT„ = a„Wi"o + ln(wiso)&nWz"o- (17) 

For the Chebyshev approximation to the flux, we know that over the entire interval v € [0, viso] the shifted Chebyshev 
polynomials have a maximum absolute value of |T'*(w)| < ±1. Therefore, the induced truncation error will be well 
approximated by the size of the coefRcient of the truncated term, i.e. 

Sc„ < A„, (18) 

due to the fact that at the LSO, the logarithmic terms are killed off. In Table I we compare the induced truncation 
errors in both the PN and Chebyshev expansions, assuming that the maximum error in both cases is at the LSO. In 
Figure ^ we present the error of the PN, Fade and Chebyshev approximations when compared to a numerical flux at 
2 to 5.5 PN order. We can see one of the downsides of the Fade method in that the 2 PN approximation has a greater 
error at the LSO than the PN flux. So in this case, the Fade approximation really only begins to make a difference 
from 3 FN onwards. We can see from this flgure that the Chebyshev approximation tries to find an equal-ripple 
error curve. It achieves this at high values of r, but as we approach the LSO the error begins to grow. Over most of 
the plotted range, the Chebyshev approximations lie below both the FN and Fade curves. In order to do this, the 
Chebyshev approximation allows the error to grow in regions where we have good agreement in order to find a better 
fit elsewhere. While the Chebyshev approximation is not as accurate as the FN or Fade at high values of r, it has no 



FIG. 1: The error of the PN, Pade and Chebyshev approximations when compared to the numerical flux for 2 to 5.5 PN order. 
The flux is modelled from R = lOOAf to the LSO at R ^ 6M. 
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FIG. 2; A logarithmic plot of the percentage error of the PN, Pade and Chebyshev approximations to the gravitational wave 
flux when compared to the numerical flux at the LSO. 



influence on our results. With the test-mass systems we have chosen, the lightest system, (100, 1A)Mq, comes into 
the detector bandwidth at about 15M, while the heaviest, (50, 1A)Mq, comes in at a greater value. We can see from 
the graph that in the region 6 < A/ < 40, in nearly all cases, the Chebyshev approximation performs better than the 
PN and the Pade at all orders of approximation. We should point out that the Chebyshev approximation is inferior 
to the Pade approximation at the 5.5 PN level. It is not understood why the Pade flux behaves so well at this order. 

We know that the main degradation in signal to noise ratio comes from a template going out of phase at the LSO. In 
Figure Q we plot the percentage errors in each of the approximations at the LSO. The PN approximation performs 
worst with errors of between 3.4% < e^^ < 42%, with the 2.5 PN approximation incurring the highest error. Besides 
the 5.5 PN order where the error is 0.06%, the error for the Pade flux at the LSO is 2.9% < e^^^^ < 7.8%. While 
we have plotted the Pade error at 5 PN order, we must remember that there is a singularity at this order, so the 
result is superfluous. In contrast, and a further sign that the Chebyshev approximation to the flux is related to the 
minimax polynomial, the error for the Chebyshev flux is approximately constant, 1.5% < e^^^^ < 1.8%, at all orders 
of approximation. This is what we would expect as it attempts to minimize the maximum error at the LSO. 



IV. RESULTS AND DISCUSSION. 



We use the technique of matched filtering to see just how well each template performs. For our 'exact' signal we 
use a restricted PN waveform where we use the exact expression for the orbital energy function, and a numerical 
gravitational wave flux function. We deflne the noise- weighted overlap between two waveforms h{t) and s{t) as the 



FIG. 3: The Fitting Factors (top) and %-error in the estimation of the chirp mass (bottom) for PN, Pade and Chebyshev 
templates as compared to an 'exact' test-mass template with parameters of (mi,m2) = (100, 1.4)Mq (left) and (mi, 7712) = 
(50, 1.4)Mq (right) . The trough in the top cell and the peak in the bottom cell for the Pade template at 5 PN is due to a 
singularity in the flux function at this order. 



inner product of the normalized waveforms denoted by 
where the scalar product is defined by 



O = , , (19) 



{h\s)^2j^^^ [hif)S*{f) + h*imf)]. (20) 

Here, an asterix denotes a complex conjugate and a tilde denotes the Fourier transform of the time domain waveform. 
Each template is a function of a number of parameters, A'^, which in turn define the dimensionality of the search 
space. For the types of systems we are considering, as well as the extrinsic parameters to and cj)o, each template 
is defined by the individual masses of the systems, (7711,7712)- Due to the short duration of the signal, parameters 
such as orbital inclination, position in the sky etc. are unimportant. While we start off with a 4-d search space, 
A'' = {^0, (f'Oj '771, "^2}, we can reduce the search to the 2-d subspace of intrinsic parameters by maximizing over the 
extrinsic parameters. Wc finally define the fitting factor J^J-^ as the overlap maximized over all parameters 



For this study, as we are concentrating on the intermediate test-mass regime, we have chosen two systems to examine 
: (100,1.4) and (50,1.4) Mq. These correspond to reduced mass ratios of 7; = 0.0136 and 0.0265 respectively. In 
order to calculate the fitting factors we used the expected EURO noise curve (as this gives us a lower frequency cutoff 
of about 10 Hz, thus giving us enough cycles to work with) |j46|. While the waveforms are described by the two 
individual mass, it is impossible for a single ground based detector to detect these parameters. What is possible is a 
combination of the individual mass called the chirp mass. This is defined as Mc — rmf'/^ . 

On the left of Figure Q we plot the fitting factors and percentage errors in the estimation of the chirp mass. The 
top cell shows the fitting factors for the PN, Pade and Chebyshev templates for the masses (100, 1.4) M©, while the 
bottom cell displays the percentage error in the chirp-mass. In general, we will take a template to be adequate if it 
reaches a predetermined threshold. This is usually a fitting factor of about 0.97 . In this case the PN templates at 
2 and 2.5 PN achieve fitting factors of about 0.96 and 0.87 respectively making them inadequate as templates. Only 
from 3 PN onwards do the templates meet the required threshold. In keeping with everything that we know about 
the PN approximation, we can see that the fitting factors follow an oscillatory pattern, i.e. the 3.5 PN is better than 
both the 3 and 4 PN templates etc. We can see from the bottom cell that not only do the lower order approximants 
achieve bad fitting factors, but they also have the highest errors in the estimation of Mc. We can also see that the 
4 PN template which performs worse than the 3.5 and 4.5 templates also has a higher parameter estimation error. 
The Pade templates achieve excellent fitting factors from 2 PN onwards, but again we can see an oscillatory nature 
in the estimation of the chirp mass. While the 2 PN template has a fitting factor of ^ 0.99 it has an error in the 
estimation of Mc almost equal to the PN template. The error in the estimation of the chirp mass begins to improve as 
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we increase the order of approximation, but we can see that the 3.5 template has a higher error in the chirp mass than 
the corresponding PN template. We should also explain that the sudden dip in the fitting factor and peak in the chirp 
mass estimation at 5 PN is due to the singularity in the flux function giving a zero fitting factor and infinite error 
in the chirp mass. At 5.5 PN order the Pade template achieves almost a perfect overlap and the lowest error in Mc- 
The Chebyshev templates, on the other hand, achieve fitting factors of almost unity at all orders of approximation. 
Not only that, but we can see that the error in estimating Mc improves, converging to a constant error value, which 
is again what we would expect as the template is near minimax. Just looking at the error in parameter estimation, 
we can see that the Chebyshev templates have achieved by 3 PN order what takes the Pade templates until 5.5 PN 
order to achieve. 

On the right hand side of Figure Q we present the same results for the masses (50, 1.4) Mq. We can see in the 
top cell that while the 2 PN template now meets the require fitting factor, the 2.5 template again falls short. Also, 
we notice that the oscillation in fitting factors is not as great as in the previous case with all PN templates achieving 
fitting factors of > 0.99 from 3 PN order onwards. Once again, the Pade templates outperform the PN templates at 
all orders. However, we again see that in order to do this, it means that the Pade templates occasionally incur larger 
errors in the estimation of the chirp mass. In this case we see that the 3 PN Pade template has an error greater than 
the corresponding PN tempalate. We also see this time that the error in the estimation of Mc is only slightly better 
in a number of cases (2 PN, 3.5 PN, 4 PN and 4.5 PN orders) than the equivalent PN templates. On the other hand, 
the Chebyshev templates again achieve fitting factors at all orders of approximation. We also find, possibly due to 
the longer waveform and the accumulation of phase information, that while the error at 2 PN order is roughly an 
order of magnitude better than the PN template, from 2.5 PN order onwards the best-fit template has parameters 
almost identical to the signal we are trying to fit. While hard to see from the plot, the error in the estimation of 
Mc is approximately 10~^% from 2.5 to 5.5 PN order. Once again, we can see that the 2.5 PN Chebyshev template 
performs as well as the 5.5 PN Pade template and better than any of the PN templates. 

V. CONCLUSION. 

In this study we have introduced a new template for detecting gravitational waves from IMRI compact binary 
systems. We have shown that one of the main problems with expansions such as the PN or Pade approximations is 
that both adhere to the Weierstrass's theorem in that we need a large number of terms to properly approximate a 
function. A better method is to use a member of the family of Ultraspherical polynomials, the Chebyshev polynomials 
of the first kind, which are closely related to the minimax polynomial due to the fact that they achieve equioscillation 
at n + 1 points in an interval. 

We demonstrated that by using the shifted Chebyshev polynomials, we can define a new Chebyshev approximation 
to the gravitational wave flux function. A major advantage of this new flux function is that when we expand a 
series in terms of Chebyshev polynomials, each subsequent coefficient is smaller than the previous one. As the shifted 
polynomials have a maximum value of unity over the domain of interest, it means that the truncation error incurred 
by going to lower orders of approximation is proportional to the neglected coefficient, which is in general a small 
number, and is much smaller than the truncation error involved in the PN approximation. 

By graphically fitting the PN, Pade and Chebyshev approximations to the gravitational wave flux function to a 
numerical flux, we were able to show that at all orders of approximation the Chebyshev approximation provided a 
better fit than the PN approximation, and was better than the Pade approximation at all orders except the 5.5 PN 
order. The closeness of the Chebyshev flux to the minimax flux was observed due to the fact that the new flux tried 
to achieve an 'equal-ripple' error curve. One of the main results of this study is the fact the new flux function has a 
lower error at the LSO, where the templates are most likely to be out of phase with a signal, than both the PN and 
Pade approximations. 

Finally, by choosing a number of test-mass systems and a fiducially 'exact' signal, we found that not only does 
the Chebyshev template always achieve higher fitting factors, but they have excellent error estimation, in most cases 
finding the parameters almost perfectly. Another of the main features of the new templates are their ability to perform 
as well at lower orders of approximation than either the PN or Pade templates at the highest order. We believe that 
these new templates will prove to be an invaluable addition to any strategy that involves the detection of gravitational 
waves using both ground and eventually, space-based detectors. 
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APPENDIX A: FITTING FACTORS AND PARAMETER ESTIMATION 

We present the values of the fitting factors and the error in the estimation of the Chirp Mass in the following tables. 
The labels T, P and C stand for the PN, Padc and Chebyshev results. Row one of each order of approximation denotes 
the fitting factor, in row two are the individual masses associated with each fitting factor and row three corresponds 
to the percentage error in the estimation of the Chirp Mass. The blank entries at 5 PN order for the Pade templates 
is due to a singularity at this order of approximation in the gravitational wave flux function. 





(100, 1.4)Mq 


(50, 1.4)M0 




T 


P C 


T P 


C 


A 

4 


0.9614 


0.9895 0.9998 


0.974 0.9863 


0.9949 




(98.83, 1.4) 


(102.13, 1.39) (99.54, 1.41) 


(49.03,1.41) (51.72,1.37) (50.08,1.4) 




0.473 


0.422 0.24 


0.37 0.083 


0.065 


tr 




0.866 


0.995 0.9995 


0.602 0.996 


0.998 




(116.19,1.3) 


(100.56, 1.4) (100.03, 1.4) 


(56.15, 1.35) (50.56, 1.39) 


(50, 1.4) 




1.62 


0.225 0.012 


2.57 0.025 





i? 
D 


0.973 


0.9969 0.9992 


0.9899 0.9985 


0.9986 




(97.52,1.41) 


(100.96, 1.39) (100.01, 1.4) 


(47.98,1.44) (50.48,1.39) (49.99,1.4) 




0.585 


0.043 0.004 


0.002 0.039 





1 


0.9975 


0.998 0.9993 


0.9955 0.9989 


0.9984 




(100.79,1.39) (100.55,1.39) (100.01,1.4) 


(50.41, 1.39) (50.21, 1.39) 


(50, 1.4) 




0.111 


0.207 0.004 


0.095 0.211 





8 


0.9843 


0.9969 0.9993 


0.9914 0.9979 


0.9984 




(103.91,1.36) (100.99,1.39) (100.01,1.4) 


(51.95,1.36) (50.49,1.39) 


(50,1.4) 




0.187 


0.031 0.004 


0.173 0.031 





9 


0.999 


0.9967 0.9993 


0.9984 0.9981 


0.9985 




(98.25, 1.42) 


(100.92, 1.39) (100.01, 1.4) 


(49.1,1.42) (50.47,1.39) 


(50, 1.4) 




0.136 


0.059 0.004 


0.107 0.047 





10 


0.996 


0.9993 


0.9957 


0.9985 




(100.88, 1.39) 


(-, -) (100.01, 1.4) 


(50.45, 1.39) (-, -) 


(50, 1.4) 




0.075 


0.004 


0.063 





11 


0.9969 


0.9999 0.9993 


0.9927 0.9995 


0.9919 




(100.78, 1.39) 


(100, 1.4) (100.01, 1.4) 


(50.43, 1.39) (50, 1.4) 


(50, 1.4) 




0.115 


0.004 0.004 


0.079 






